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Abstract 

A  passive  liquid-feed  direct  methanol  fuel  cell  (DMFC)  with  neither  external  liquid  pumps  nor  gas  blowers  is  modeled  mathematically. 
A  one-dimensional  model  is  developed  by  considering  inherently  coupled  heat  and  mass  transport,  along  with  the  electrochemical  reactions 
occurring  in  passive  DMFCs.  The  analytical  solutions  predicting  the  performance  of  this  type  of  fuel  cell  operating  with  different  methanol 
concentrations  are  obtained.  It  is  shown  that  the  performance  of  passive  DMFCs  increases  with  methanol  concentration.  It  is  further  revealed 
that  the  improved  performance  with  higher  methanol  concentrations  is  due  primarily  to  the  increased  operating  temperature  resulted  from  the 
exothermic  reaction  between  permeated  methanol  and  oxygen  on  the  cathode. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  an  electro¬ 
chemical  device  that  directly  converts  chemical  energy  into 
electrical  energy.  Because  of  its  high  efficiency,  high  energy 
density,  low  emission,  and  simple  structure,  the  DMFC  has 
been  projected  to  be  a  prime  candidate  for  powering  portable 
devices  [1].  In  terms  of  fuel  and  oxidant  delivery  schemes, 
DMFCs  can  be  classified  either  as  active-feed  DMFCs  or  as 
passive-feed  DMFCs. 

In  active-feed  DMFCs,  methanol  solution  is  usually 
delivered  by  a  liquid  pump,  while  oxidant  (oxygen  or  air)  is 
supplied  by  a  gas  blower/fan.  To  improve  cell  performance, 
most  previous  studies  on  active-feed  DMFCs  have  been 
focused  on  the  following  aspects:  improving  the  electro¬ 
activity  of  methanol  oxidation  on  the  anode  by  exploring 
more  active  electrocatalysts  [2-4];  optimizing  cathode 
electrode  structures  to  avoid  severe  flooding  [5];  investi¬ 
gating  the  effect  of  methanol  crossover  [6-11];  reducing 
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the  substantial  methanol  crossover  by  modifying  existing 
polymer  membranes  [12]  or  searching  for  alternatives  [13]; 
flow  field  designs  and  CO2  management  [14-19]. 

In  contrast,  the  concept  of  passive  DMFCs  is  that  a  DMFC 
can  operate  stand-alone,  eliminating  the  external  means  of 
liquid  transport  and  air  movement.  Therefore,  the  fuel  cell 
system  becomes  much  simpler  and  more  compact.  This  con¬ 
cept  provides  a  unique  feature  for  the  miniaturization  of  fuel 
cells.  For  this  reason,  passive  DMFCs  have  recently  received 
much  attention.  Kim  et  al.  [20]  fabricated  and  tested  a  sin¬ 
gle  cell  and  monopolar  DMFC  stack  operating  under  passive 
and  air-breathing  conditions.  Liu  et  al.  [21]  studied  sintered 
stainless  steel  fiber  felt  as  the  gas  diffusion  layer  in  an  air- 
breathing  DMFC.  Shimizu  and  Momma  [22]  reported  their 
activities  regarding  the  research  and  development  of  DMFCs 
that  operated  passively  at  room  temperature.  Chen  and  Yang 
[23]  investigated  the  effect  of  operating  conditions  on  the 
power  density  of  an  air-breathing  DMFC. 

As  compared  with  their  active  counterparts,  passive- 
feed  DMFCs  have  much  lower  power  density  (about 
20mWcm-2)  because  of  two  major  problems.  One  is  the 
inability  to  handle  the  excess  water  evolved  on  the  cathode 
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and  crossed  from  the  anode.  The  other  is  that  conventional 
designs  lose  too  much  heat  from  the  fuel  cell  to  the  ambient 
air,  resulting  in  rather  low  cell  operating  temperatures  (less 
than  30  °C).  Therefore,  the  ambient  conditions  under  which 
passive  DMFCs  operate  require  added  attention  to  be  given 
to  thermal  and  water  management. 

Because  several  interrelated  processes  occur  simultane¬ 
ously,  the  development  of  a  mathematical  model  is  critical  to 
the  design  and  the  optimization  of  passive-feed  DMFCs.  In 
particular,  understanding  thermal  effects,  which  are  insepa¬ 
rable  from  the  management  of  water  in  the  cell,  is  necessary 
to  optimize  the  performance  of  these  fuel  cells.  A  number 
of  studies  [24-34]  have  been  reported  simulating  PEM  fuel 
cells  and  DMFCs,  but  most  of  them  were  isothermal  mod¬ 
els,  except  for  a  few  papers  that  took  into  account  thermal 
effects  for  simulating  solid  oxide  fuel  cells  [35-37].  Nord- 
lund  and  Lindbergh  [30]  proposed  an  isothermal  agglomerate 
model  based  on  the  reaction  mechanism  for  the  electrochemi¬ 
cal  oxidation  of  methanol  to  study  the  influence  of  the  porous 
structure  on  the  direct  methanol  fuel  cells.  Wang  and  Wang 
[31]  used  two-phase  multi-component  model  to  simulate  a 
DMFC.  The  anode  and  cathode  electrochemical  reactions, 
diffusion  and  convection  of  both  gas  and  liquid  phases  in 
the  backing  layers  and  flow  channels,  mixed  potential  effect 
due  to  methanol  crossover  and  the  effect  of  methanol  feed 
concentration  were  explored.  Murgia  et  al.  [32]  developed  a 
one-dimensional,  two-phase,  multi-component  steady-state 
model  based  on  the  phenomenological  transport  equations. 
In  their  model,  the  influences  of  the  methanol  concentration, 
the  pressure  gradient  in  the  catalyst  layer,  flooding  in  the  cath¬ 
ode  gas  diffusion  layer,  and  methanol  crossover  on  the  cell 
performance  were  investigated. 

In  this  paper,  we  present  a  theoretical  model  that  incorpo¬ 
rates  the  effects  of  coupled  heat  and  mass  transfer  and  electro¬ 
chemical  kinetics  in  a  passive-feed  DMFC.  We  demonstrate 
that  this  simplified  model  is  capable  of  shedding  light  on  the 
mechanisms  leading  to  a  better  performance  when  this  type  of 
fuel  cells  operate  with  higher  methanol  concentrations  under 
ambient  conditions. 


2.  Model  development 

Consider  a  passive-feed  DMFC  shown  in  Fig.  1,  which 
consists  of  a  fuel  tank,  an  anode  current  collector  (ACC), 
an  anode  gas  diffusion  layer  (AGDL),  an  anode  catalyst 
layer  (ACL),  a  membrane,  a  cathode  catalyst  layer  (CCL),  a 
cathode  gas  diffusion  layer  (CGDL),  and  a  cathode  current 
collector  (CCC).  On  the  anode,  diluted  methanol  solution 
is  introduced  to  the  reaction  zone  without  any  external 
means  of  liquid  transport.  From  the  ACC  through  the  AGDL 
and  finally  to  the  ACL,  methanol  solution  is  transported 
primarily  by  diffusion.  After  the  electrochemical  reaction 
of  methanol  oxidation,  CO2  produced  at  the  ACL  moves 
counter-currently  toward  the  fuel  tank  via  the  AGDL  and 
ACC.  At  sufficiently  high  current  densities,  CO2  emits  in 


Fig.  1.  Schematic  of  a  passive  liquid-feed  DMFC. 


the  form  of  gas  bubbles  from  the  surface  of  the  ACC.  For 
the  vertically  orientated  DMFC  shown  in  Fig.  1,  in  the  fuel 
tank,  due  to  buoyancy,  CO2  moves  upwards  adjacent  to  the 
surface  of  the  ACC,  while  the  liquid  methanol  solution  away 
from  the  surface  of  the  ACC  moves  downwards.  As  such,  the 
diffusion  of  methanol  solution  from  the  fuel  tank  to  the  ACC 
is  virtually  enhanced  by  natural  convection.  On  the  cathode 
of  this  passive-feed  DMFC,  oxygen  is  taken  passively  from 
the  ambient  air  without  any  means  of  air  movement.  Almost 
all  the  heat  generated  by  the  fuel  cell  is  lost  from  the  vertical 
CCC  to  the  ambient,  leading  to  natural  heat  convection  at 
the  surface  of  the  CCC.  The  mass  transport  of  oxygen  and 
water  on  the  cathode  is  enhanced  by  natural  convection. 

It  is  clear  from  the  above  description  that  the  operation 
of  the  passive-feed  DMFC  considered  in  this  work  actually 
involves  rather  complex  transport  processes  of  heat  and  mass, 
which  are  linked  to  the  electrochemical  reactions  of  methanol 
oxidation  on  the  anode  and  oxygen  reduction  on  the  cathode. 
To  make  the  complicated  process  tractable,  we  treat  it  as  a 
one-dimensional  problem,  with  the  v-axis  origin  set  at  the 
ACL,  as  illustrated  in  Fig.  1.  We  also  introduce  the  following 
simplifications  and  assumptions: 

(1)  The  fuel  cell  is  assumed  to  operate  under  steady-state 
conditions. 

(2)  Compared  with  the  heat  generated  by  electrochemical  re¬ 
action  and  overpotential,  Joule  heat  caused  by  the  current 
flow  through  each  component  is  ignored. 

(3)  The  fuel  tank,  ACC,  and  AGDL  are  well  insulated  from 
the  ambient.  No  heat  is  lost  from  these  components  to 
the  ambient.  Thus,  the  temperature  at  each  of  these  com¬ 
ponents  is  the  same  as  that  at  ACL. 

(4)  Considering  that  ACL  and  CCL  are  much  thinner  than 
AGDL,  CGDL  and  the  membrane,  they  are  treated  as  an 
interface  and  the  temperature  and  concentration  distribu¬ 
tions  are  uniform. 

(5)  The  transport  of  heat  and  mass  through  the  gas  diffusion 
layers  is  assumed  to  be  a  diffusion-predominated  process 
and  the  convection  effect  due  to  bulk  flow  is  ignored. 
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(6)  On  the  cathode,  the  heat  and  mass  are  transferred  between 
the  CCC  and  the  ambient  by  natural  convection. 

(7)  On  the  anode,  the  mass  transfer  of  methanol  from  the 
bulk  solution  to  the  ACC  is  assumed  to  be  driven  by 
natural  convection.  For  a  sufficiently  large  fuel  tank,  the 
methanol  concentration  in  the  tank  during  the  fuel  cell 
discharging  is  assumed  to  be  constant. 

(8)  Methanol  crossover  through  the  membrane  is  assumed 
to  be  due  to  the  combined  effect  of  the  concentration 
gradient  between  the  anode  and  cathode  and  the  electro¬ 
osmosis  force. 

(9)  Since  the  reaction  rate  of  methanol  on  the  cathode  is 
rather  fast,  the  methanol  concentration  in  the  CCL  is  as¬ 
sumed  to  be  depleted. 

2.7.  Mass  transport 


where  and  Cj£f  represent,  respectively,  the  reference  ex¬ 
change  current  density  on  the  anode  and  the  reference  concen¬ 
tration  of  methanol,  and  Cm,aci  is  the  methanol  concentration 
in  the  ACL.  The  permeation  flux  of  methanol  through  the 
membrane,  AcroSs,  can  be  determined  from: 

AT  —  —  Dcff  d^m,mem  .  m  J_  n  , 

iVcross  ^m,mem  ^  '"'dp’  G/ 

where  2)^fmem  is  the  effective  diffusion  coefficient  of 
methanol  in  the  membrane  and  n ™  is  the  electro-osmotic  drag 
coefficient  of  methanol. 

On  the  cathode,  oxygen  extracted  from  the  air  reacts  with 
the  electron  and  proton  to  produce  water.  Additionally,  part 
of  oxygen  is  consumed  due  to  methanol  crossover  to  form  an 
internal  current  and  a  mixed  potential.  Therefore,  the  oxygen 
flux  from  the  CGDL  to  the  CCL  can  be  determined  from: 


The  transport  process  of  methanol  from  the  fuel  tank  to 
the  ACC  is  described  by: 


Nm  —  ^m(Gm5tank  acc^’ 


m,aco 


(1) 


where  Nm  represents  the  methanol  flux,  hm  the  mass  transfer 
coefficient  at  the  ACC  surface,  and  acc  and  Cm?tank  are  the 
methanol  concentrations  at  the  ACC  surface  and  in  the  tank, 
respectively. 

In  the  ACC  and  AGDL,  since  no  electrochemical  reaction 
occurs,  the  methanol  flux  remains  the  same  and  is  related  to 
the  concentration  gradient  by  Fick’s  law  as: 


Nm  =  ~D 


eff 
m,acc 


dC 


m,acc 


dx 


and 


Nm  =  ~D 


eff 


dC 


m,agdl 


m'a§dl  dx 


where  D^facc  and  D^agdl  represent  the  effective  diffusion 
coefficients  of  methanol  in  the  ACC  and  AGDL,  respectively. 
Combining  Eqs.  (1)— (3)  yields: 


1  .  3 

TTGc  +  “  Across 

4  F  2 


To  account  for  the  effect  of  methanol  crossover  on  the  cathode 
overpotential,  we  assume  that  the  methanol  from  the  anode 
completely  reacts  electrochemically  on  the  cathode.  Based  on 
this  assumption,  the  internal  current,  /p,  due  to  the  methanol 
oxidation  on  the  cathode  can  be  obtained  based  on  the  per¬ 
meation  flux  of  methanol  [31]: 


—  6FN cross  • 


It  follows  from  Eqs.  (6)  and  (9)  that  the  cathode  overpotential 
taking  account  of  the  mixed  potential  can  be  determined  from: 


(10) 


where  and  represent,  respectively,  the  reference  ex¬ 
change  current  density  on  the  cathode  and  the  reference  con¬ 
centration  of  oxygen,  and  Co2,ccl  is  the  oxygen  concentration 
in  the  CCL. 

In  the  CGDL  and  CCC,  the  oxygen  flux  remains  the  same 
and  is  related  to  the  concentration  gradient  by: 


Nm  —  Ot\(C  m,t  ank  ^m,SLC  l), 


(4) 


-l 


where  ot\  = 


*agdl 


D 


eff 


+ 


!acc 


D 


eff 


+ 


h 


m,agdl  ^m,acc 

The  methanol  flux,  Nm,  is  related  to  the  proton  current 
density,  za,  and  the  permeation  flux  of  methanol  through  the 
membrane,  Ncross,  by: 


1 

N in  —  T  N, 


6  F 


cross  > 


(5) 


where  F  is  the  Faraday’s  constant  (96,485  C  mol-1).  The 
current  density,  /a,  can  be  determined  from  the  Tafel 
equation: 


i 


a 


.m  Wn,ac  1 
Lef  ~Fm 
We  f 


exp 


/  a^F 
\  7?raci 


No2  = 
and 


—  D 


eff  dCo2?cgdl 
°2’cgdl  dx 


No,  =  ~D 


eff 
C>2,ccc 


dCo2, 


CCC 


dx 


(11) 


(12) 


where  77 q1  ccrdl  and  Dq1  ccc  represent,  respectively,  the  effec¬ 
tive  diffusion  coefficients  of  oxygen  in  the  CGDL  and  in  the 
CCC. 

The  oxygen  transport  from  the  ambient  air  to  the  CCC  is 
due  to  natural  convection  and  can  be  described  by: 

Nq2  =  ho2  (Cq2,  amb  —  ^C>2,ccc)’  (13) 


where  ho2  is  the  mass  transfer  coefficient  of  oxygen  on  the 
cathode,  Cqo  ccc  and  Cq2,  amb  are  the  oxygen  concentrations 
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at  the  surface  of  the  CCC  and  in  the  bulk  air,  respectively. 
Combining  Eqs.  (1 1)— (13)  gives  the  oxygen  flux  as: 


Nq2  —  <^2(Co2,amb  Cq2,  ccl)> 


(14) 


where  a  2  = 


+ 


*cgdl 


D 


,eff 

02,cgdl 


+ 


^ccc 


D 


eff 

02,ccc 


-1 


2.2.  Heat  transport 


On  the  anode,  heat,  generated  by  the  electrochemical  re¬ 
action  in  the  ACL,  is  given  by: 


^acl  —  l 


A//a- A Ga\ 
nF  )  ’ 


(15) 


accumulated  on  the  cathode  due  to  the  limited  removal  rate 
of  water  with  no  external  means  of  air  movement.  Then,  the 
flux  of  water  evaporation  can  also  be  determined  by  natural 
convection: 

Nh20  =  ^H2o(^h20,ccc  -  ^H20,amb),  (20) 

where  /*h2o  denotes  the  mass  transfer  coefficient  of  water 
vapor  on  the  cathode,  ccc  and  CH2o,amb  are  the  water 
vapor  concentrations  at  the  surface  of  the  CCC  and  in  the  bulk 
air,  respectively.  The  water  vapor  concentration  can  be  deter¬ 
mined  fromPH2o/^7\  and  the  saturated  pressure  of  moist  air 
is  determined  from  [27]: 

d  _  1  n(-2. 1794+0.02953  x  T-9. 1837  x  1(T5  x  T2  +  1.4454x  10“7  x  T3) 

^sat — AU 


where  the  first  term  represents  the  heat  due  to  the  activation 
and  mass  transfer  overpotentials  on  the  anode,  the  second 
term  represents  the  entropy  change  of  the  anodic  electro¬ 
chemical  reaction,  with  A//a  denoting  the  anodic  reaction 
enthalpy  and  AGa  the  Gibbs  free  energy.  Eq.  (15)  can  be 
rewritten  as  [38]: 


4aci  =  iri  a  -Pi-  ft(^a  cl  -  298),  (16) 

where  ft  =  ^=(A H®  -  AG°)  and  =  pp(cco2  -  Cm  - 
Ch2o)- 

Neglecting  Joule  heat  in  the  membrane,  we  can  relate  the 
heat  flux,  gac i,  to  the  temperature  gradient  across  the  mem¬ 
brane  as: 


dT 

<7acl  —  2,mem—  , 

dx 


(17) 


where  kmem  is  the  effective  thermal  conductivity  of  the  mem¬ 
brane. 

The  heat  generation  on  the  CCL,  qcc i,  can  be  determined 
from: 


A Hc  -  AGC 

<7ccl  =  0  +  lv)0c  -  I - - - ftAft2o,  (18) 

where  the  first  term  represents  the  heat  due  to  activation  and 
mass  transfer  overpotentials  and  mixed  potential  caused  by 
methanol  crossover  on  the  cathode;  the  second  term  accounts 
for  the  entropic  loss,  with  A Hc  denoting  the  cathodic  reaction 
enthalpy  and  A Gc  is  the  Gibbs  free  energy;  the  third  term 
reflects  the  heat  due  to  the  evaporation  of  liquid  water  at  the 
CCL,  with  hv  denoting  the  latent  heat  of  water  and  Nu2o  is 
the  flux  of  water  evaporation.  Eq.  (18)  can  be  rewritten  as: 


qcc\  =  (i  +  iP)ric  -  ft  -  ft(ft cl  -  298)  -  hwNu2o ,  (19) 

where  £3  =  37 (A H®  -  A G°)  and  Pa  =  jp  (Ch2o- 

lCo2). 

We  assume  that  the  air  at  the  CCC  is  in  a  saturated  state 
based  on  the  following  justifications:  in  addition  to  the  water 
continuously  generated  on  the  cathode,  the  permeation  rate 
of  water  from  the  anode  to  cathode  is  relatively  high  in  DM- 
FCs,  even  under  the  open  circuit  voltage  (OCV)  condition 
[34].  In  particular,  in  passive  DMFCs,  more  water  may  be 


(21) 


The  total  heat  generation  qiot  is  given  by: 

qtot  =  Qdid  T  ^ccl? 


(22) 


which  is  related  to  the  temperature  gradient  in  the  CGDL  and 
the  CCC,  respectively,  to 


^tot  —  ftgdl 


d  T 

dx  ’ 


(23) 


and 


qtot 


=  -x 


ccc 


dr 

dx  ’ 


(24) 


where  Acgcu  and  Accc  represent  the  effective  thermal  conduc¬ 
tivities  of  in  the  CDGL  and  in  the  CCC.  Finally,  the  heat 
transfer  from  the  CCC  to  the  ambient  air  can  be  described 
using  the  Newton’s  cooling  law  as: 

qtot  —  ht(Tccc  —  ramb),  (25) 


where  h t  is  the  heat  transfer  coefficient  due  to  natural  con¬ 
vection,  Tc°cc  and  Tamb  are  the  temperatures  at  the  surface  of 
the  CCC  and  the  ambient  air.  Combining  Eqs.  (23)-(25),  we 
obtain: 


^tot  —  oi^(Tcc\  ramb),  (26) 

where  <23  =  ( j-  + 

\  h[  ^-cgdl  Acc  / 

The  heat  transfer  coefficient  at  the  surface  of  the  CCC  can 
be  determined  from  [39]: 


hL  0.67Ra\/4 

Nu  =  —  =  0.68  + - L 

A  [1  +  (0.492/ PrC16]/9 


(27) 


where  Ra l  is  the  Rayleigh  number  ( RaL  =  GrPr ),  Pr 
the  Prandtl  number,  and  Gr  is  the  Grashof  number 
^Gr  =  fjg(Tw-To)L  y  Based  on  the  analogy  principle  be¬ 
tween  heat  and  mass  transfer,  the  mass  transfer  coefficient 
can  also  be  determined  from  Eq.  (27)  with  Nu  replaced  by 
Sherwood  number,  Sh  =  /?mgsL ,  and  Pr  replaced  by  Lewis 
number,  Le  =  Sh/Nu  (for  gas  Le=  1  [40]).  On  the  anode,  the 
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natural  convection  is  caused  by  bubble  behavior.  The  mass  where  £%  is  the  open  circuit  voltage  at  T=  298  K  and 

transfer  coefficient  at  the  ACC  is  determined  by  assuming  (dE/dT)  represents  the  rate  of  change  of  electromotive  force, 

that  Sh/Nu  =  2.0. 


2.3.  Cell  performance 


3.  Analytical  solutions 


With  the  methanol/oxygen  concentration  at  the  catalyst 
layers  and  the  temperature  distributions,  and  the  anodic  and 
cathodic  overpotentials  obtained  from  the  equations  pre¬ 
sented  above,  we  can  assess  the  cell  performance  based  on: 

Well  —  ^cell  Osl  0c  ^celh  (28) 


We  now  obtain  analytical  solutions  from  the  equations 
presented  in  the  preceding  section.  Combining  Eqs.  (4),  (5), 
and  (7)  yields  the  methanol  concentration  at  the  ACL: 
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(30) 


where  Lceii  is  the  thermodynamic  equilibrium  potential  of  the 
fuel  cell  and  is  a  function  of  temperature  and  pressure,  and 
Rcq\\  is  the  internal  resistance  of  the  fuel  cell.  The  thermody¬ 
namic  equilibrium  potentials  of  the  fuel  cell  can  be  calculated 
by: 

Seen  =  £c°ell  +  AT  (!jf )  .  (29) 


With  the  aide  of  Eq.  (30),  the  oxygen  concentration  at  the 
CCL  can  be  determined  by  using  Eqs.  (8)  and  (14)  to  give: 
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Table  1 


Physical  properties 


Parameter/symbol  (unit) 

Value 

Refs. 

Liquid  methanol  enthalpy  of  formation///m  (Jmol-1) 

-238.66  x  103 

Liquid  water  enthalpy  of  formation/ AHW  (J  mol-1) 

-285.83  x  103 

Carbon  dioxide  enthalpy  of  formation/ A  Hco2  (J  mol-1) 

-393.51  x  103 

Liquid  methanol  Gibbs  free  energy/ A  Gm  (J  mol-1) 

-166.27  x  103 

Liquid  water  Gibbs  free  energy/ A  Gw  (J  mol-1) 

-237.08  x  103 

Carbon  dioxide  Gibbs  free  energy/ A  Gqo2  (J  mol-1) 

-394  x  103 

Liquid  methanol  specific  heat  capacity /cm  (J  mol-1  K-1) 

80.96 

Liquid  water  specific  heat  capacity /cw  (J  mol- 1  K- 1 ) 

75.24 

Carbon  dioxide  specific  heat  capacity/cco2  (J  mol-1  K-1) 

36.9 

Oxygen  specific  heat  capacity/ccb  (J  mol-1  K-1) 

39.44 

Evaporation  heat  of  water thy  (J  mol-1) 

44.86  x  103 

Pressure  of  air  in  cathode/Po  (Pa) 

1  x  105 

Oxygen  concentration  in  the  ambient  air/Co2,amb  (molm-3) 

0.21  x  PIRT 

Thickness  of  ACC  and  CCC//acc  (m) 

0.001 

Assumed 

Thickness  of  AGDL  and  CGDL//agdi  (m) 

0.0003 

[41] 

Electro-osmotic  drag  coefficient  of  water/nHoO 

2.9exp(1029/(l/333  —  1/7)) 

[42] 

Electro-osmotic  drag  coefficient  of  methnaol/n™ 

xm  x  «h2  0 

[31] 

Conductivity  of  membraneArm  (S  m- 1 ) 

7.3  exp(1268(l/298  —  l/T)) 

[43] 

Correction  factor  of  diffusion  coefficient^: 

0.6 

Assumed 

Thermal  conductivity  of  mernbrane/A.mcm  (W  m-1  K-1) 

0.21 

[33] 

Thermal  conductivity  of  CGDL/kcgdi  (W  m-1  K-1) 

1.6 

[44] 

Thermal  conductivity  of  CCC/A.ccc  (W  m- 1  K- 1 ) 

16 

Diffusion  coefficient  of  methanol  in  mernbrane/D^fmem  (m2  s-1) 

4.9  x  10“ 10  exp(2436(l/333  -  1/7)) 

[43] 

Diffusion  coefficient  of  methanol  in  ACC/D^facc  (m2  s_1) 

2.8  x  10“9 

[34] 

Diffusion  coefficient  of  methanol  in  AGDL/D^fagdl  (m2  s-1) 

2.8  x  10“9£ 

Diffusion  coefficient  of  oxygen  in  CCC  and  air/Z)QfJfccc  (m2  s-1) 

0.25  x  10"4 

Diffusion  coefficient  of  oxygen  in  CGDL/Dq^  cgdl  (m2  s-1) 

0.25  x  10-4s 

Diffusion  coefficient  of  water  vapor  in  air/D^0  air  (m2  s-1) 

6  x  10“4 

Assumed 

Transfer  coefficient  of  anode/aa 

0.35 

Assumed 

Transfer  coefficient  of  cathode/ac 

0.8 

Assumed 

Reference  concentration  of  methanol/C™f  (molm-3) 

4  x  103 

Assumed 

Reference  concentration  of  oxygen/C^  (mol/m3) 

0.21  x  P()/RT 

Reference  exchange  current  density  on  anode/iJJf  (Am-2) 

94.25  exp(35,570/P(l/353  -  l/T)) 

[31] 

Reference  exchange  current  density  on  cathode//^  (Am-2) 

0.04222  exp(73,200/P(l/353  -  l/T)) 

[31] 

Rate  of  change  of  electromotive  force/3E/3T  (VK-1) 

-1.4  x  10“4 

[34] 
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Utilizing  the  calculated  methanol  and  oxygen  concentrations 
at  the  catalyst  layer,  the  temperatures  at  catalyst  layer  are 
obtained  from  Eqs.  (16),  (17),  (19),  (22),  and  (26).  The  tem¬ 
peratures  in  the  ACL  and  the  CCL  is  given,  respectively,  as: 


variation  in  the  temperature  at  the  ACL  with  current  density 
for  different  methanol  concentrations.  The  increased  oper¬ 
ating  temperature  improves  the  electrochemical  reactions  of 


if]  a  ~  Pi  +  29802  +  i?c(i  +  ip)  ~P3  +  298 04  +  f^(ir]apA  -  Pi  Pa  +  2980204  -  «30i  +  29802a3 

r  Amem 

T  Tamb  —  hyNn2o 

«3  +  Pi  +  PA  +  fSSL(p20l3  +  0204) 

Amem 


(32) 


and 

Tccl  =  ^ad  -  ~~—~(il1a  -Pl~  Tael 02  +  298 02).  (33) 

^mem 

It  should  be  mentioned  that  the  heat  and  mass  transfer  coef¬ 
ficients  are  determined  based  on  the  iteration  procedure  by 
presuming  a  wall  temperature  at  the  CCC  surface.  Iteration 
will  not  cease  until  a  right  wall  temperature  and  right  heat, 
and  mass  transfer  coefficients,  obtained  from  Eq.  (27),  are 
determined  such  that  the  total  heat  generated  is  balanced  by 
the  heat  loss  from  the  cathode  to  the  ambient.  The  results  to 
be  presented  in  the  next  section  were  obtained  based  on  the 
parameters  listed  in  Table  1 . 


4.  Results  and  discussion 

Lig.  2  shows  the  temperature  distribution  in  the  cell  when 
it  is  filled  with  2  M  methanol  solution  and  operates  at  different 
current  densities.  The  four  data  points  represent,  respectively, 
the  temperatures  at  the  ACL,  the  CCL,  the  CGDL/CCC  inter¬ 
face,  and  the  outer  surface  of  the  CCC.  It  is  seen  from  Lig.  2a 
that  at  a  lower  current  density  the  temperature  on  the  anode  is 
lower  than  that  on  the  cathode.  This  is  because  the  heat  gener¬ 
ation  rate  by  the  anodic  overpotential  at  low  current  densities 
is  less  than  the  endothermic  heat  demanded  by  the  electro¬ 
chemical  reaction  of  methanol  oxidation.  As  a  result,  some 
heat  has  to  be  taken  from  the  cathode.  With  an  increase  in 
current  density,  the  heat  generation  rate  by  the  anodic  over¬ 
potential  increases.  Therefore,  as  indicated  in  Lig.  2b,  the 
temperature  on  the  anode  becomes  higher  than  that  on  the 
cathode.  With  a  further  increase  in  current  density,  the  tem¬ 
perature  difference  between  the  anode  and  cathode  increases, 
as  evident  by  comparing  Lig.  2b  and  c. 

Lig.  3  shows  the  polarization  curves  when  the  fuel  cell 
operates  with  1  M,  2  M,  3  M,  and  4  M  methanol  solutions,  re¬ 
spectively.  It  can  be  seen  from  this  figure  that  with  an  increase 
in  methanol  concentration,  the  cell  performance  upgrades 
progressively.  A  higher  methanol  concentration  results  in  a 
higher  methanol  permeation  rate  from  the  anode  side  to  the 
cathode  side.  On  the  cathode  side,  methanol  reacts  with  the 
oxygen  electrochemically  to  form  a  mixed  potential.  Hence, 
a  higher  methanol  concentration  leads  to  a  higher  mixed  po¬ 
tential,  thereby  causing  a  higher  heat  generation  rate.  As  a 
result,  a  higher  concentration  results  in  a  higher  operating 
temperature.  This  is  evident  from  Lig.  4,  which  shows  the 


both  methanol  oxidation  and  oxygen  reduction,  producing 
a  higher  voltage  for  a  given  current  density.  Therefore,  it 
can  be  concluded  that  the  improved  cell  performance  with 
a  higher  methanol  concentration  is  primarily  attributed  to 


x  (m) 

(c)  Current  density  i=50mA/cm2 


Fig.  2.  Temperature  distribution  in  the  cell  at  different  current  densities:  (a) 
current  density,  i  -  1  mA  cm-2,  (b)  i  =  5  mA  cm-2,  and  (c)  i  =  50  mA  cm~2. 
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Fig.  3.  Polarization  curves  of  the  passive  DMFC  operating  with  different 
methanol  concentrations. 

the  increased  operating  temperature  caused  by  the  exother¬ 
mic  reaction  between  permeated  methanol  and  oxygen  on  the 
cathode.  It  is  also  noted  from  Fig.  4  that  for  each  methanol 
concentration,  the  temperature  at  the  ACL  increases  with  cur¬ 
rent  density.  This  is  because  the  heat  generation  rate  due  to 
activation  and  mass  transfer  activation  increase  with  an  in¬ 
crease  in  current  density. 

In  order  to  further  verify  that  the  improved  cell  perfor¬ 
mance  with  higher  methanol  concentrations  is  due  to  the  in¬ 
creased  operating  temperature  caused  by  higher  methanol 
permeation  rates,  we  now  let  the  passive  DMFC  still  oper¬ 
ate  with  4  M  methanol,  but  remove  the  excess  heat  generated 
from  the  system  by  cooling  down  the  cell  to  a  lower  tem¬ 
perature.  Fig.  5  presents  the  performance,  represented  by  the 
dashed  curve,  of  the  fuel  cell  operating  with  4  M  methanol 
and  at  an  operating  temperature  controlled  to  be  the  same 
as  that  when  the  cell  operates  with  2  M  methanol  under  the 
open  circuit  voltage  condition.  For  comparison,  the  polar¬ 
ization  curves  with  2M  and  4M  without  cooling  are  also 
shown  in  Fig.  5.  It  can  be  seen  from  this  figure  that  feed¬ 
ing  4  M  methanol  solution  with  a  lowed  operating  tempera¬ 
ture  results  in  significantly  lower  performance  than  feeding 


Fig.  4.  Variations  in  the  operating  temperature  (at  the  anode  catalyst  layer) 
with  current  density  for  different  methanol  concentrations. 


Fig.  5.  Polarization  curves  for  the  4  M  methanol  operation  with  and  without 
the  operating  temperature  lowered  to  the  value  of  the  2  M  methanol  operation 
under  the  OCV  condition. 

the  same  concentration  methanol  without  cooling  does.  Al¬ 
though  the  performance  with  4  M  methanol  at  the  lowered 
operating  temperature  becomes  much  lower  than  that  with¬ 
out  the  control  of  the  operating  temperature,  it  is  still  a  bit 
higher  than  that  with  2  M  methanol.  This  may  be  due  to  the 
higher  mass  transfer  rate  with  a  higher  methanol  concentra¬ 
tion.  The  results  presented  in  Fig.  5  indicate  that  the  excess 
heat  generated  by  the  exothermic  reaction  of  the  permeated 
methanol  on  the  cathode  is  removed  from  the  cell,  the  cell  per¬ 
formance  declines  and  becomes  almost  the  same  as  that  with 
2  M  methanol.  The  comparison  in  cell  performance  among 
the  three  operating  conditions  confirms  that  the  improved  per¬ 
formance  of  a  passive  DMFC  running  with  higher  methanol 
concentrations  is  attributed  primarily  to  the  higher  operating 
temperature  caused  by  the  exothermic  reaction  between  the 
permeated  methanol  and  oxygen  on  the  cathode. 

The  effect  of  room  temperature  on  the  cell  performance 
with  different  methanol  concentrations  is  shown  in  Fig.  6.  It 
is  clear  from  this  figure  that  the  cell  performance  increases 
with  room  temperature.  The  lower  heat  loss  rate  to  the  am¬ 
bient  with  a  higher  room  temperature  creates  a  higher  oper¬ 
ating  temperature,  causing  the  electrochemical  kinetics  of 


Fig.  6.  Polarization  curves  for  different  methanol  concentrations  and  at  dif¬ 
ferent  room  temperatures. 
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both  methanol  oxidation  and  oxygen  reduction  to  be  im¬ 
proved.  The  results  presented  in  Fig.  6  suggest  that  the  perfor¬ 
mance  of  a  passive  DMFC  is  sensitive  to  a  change  in  ambient 
conditions. 

5.  Concluding  remarks 

A  one-dimensional  thermal  model  for  passive  DMFCs  is 
developed.  Inherently  coupled  heat  and  mass  transport,  along 
with  the  electrochemical  reactions,  are  considered.  The  ana¬ 
lytical  solutions  of  methanol  and  oxygen  concentration,  and 
temperature  distributions  through  the  cell  components  are 
obtained.  The  performance  of  the  fuel  cell  fed  with  different 
methanol  concentrations  (1-4  M)  is  predicted.  It  has  been 
shown  that  at  low  current  densities,  the  temperature  on  the 
anode  is  lower  than  that  on  the  cathode.  With  an  increase  in 
current  density,  the  temperature  on  the  anode  becomes  higher 
than  that  on  the  cathode.  A  higher  current  density  leads  to  a 
higher  temperature  difference  between  the  anode  and  cath¬ 
ode.  Furthermore,  with  this  model  it  has  been  revealed  that 
passive  DMFCs  operating  with  a  higher  methanol  concentra¬ 
tion  yield  a  better  performance.  The  improved  performance 
with  higher  methanol  concentrations  has  been  found  to  be 
due  primarily  to  the  higher  operating  temperature  caused  by 
the  exothermic  reaction  between  the  permeated  methanol  and 
oxygen  on  the  cathode.  The  results  presented  in  this  work  sug¬ 
gest  that  thermal  management  in  passive  DMFCs  is  critical 
for  improving  the  performance  of  this  type  of  fuel  cells. 
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